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Dewetting of Thin Viscoelastic Polymer Films on Slippery Substrates 
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Dewetting of thin polystyrene films deposited onto silicone wafers at temperatures close to the 
glass transition exhibits unusual dynamics and front morphologies. Here, we present a new theoreti- 
cal approach of these phenomena taking into account both the viscoelastic properties of the film and 
the non-zero velocity of the film at the interface with the substrate (due to slippage) . We then show 
how these two ingredients lead to : (a) A very asymmetric shape of the rim as the film dewetts, 
(b) A decrease of the dewetting velocity with time like t" for times shorter than the reptation 
time (for larger times, the dewetting velocity reaches a constant value). Very recent experiments 
by Damman, Baudelet and Reiter [Phys. Rev. Lett. 91, 216101 (2003)] present, however, a much 
faster decrease of the dewetting velocity. We then show how this striking result can be explained 
by the presence of residual stresses in the film. 

PACS numbers: 68.60.-p, 68.15.+e, 68.55.-a, 83.10.-y 



Thin liquid films are of great scientific and technolog- 
ical importance, and display a variety of interesting dy- 
namics phenomena 0,0- In engineering, for instance, 
they serve to protect surfaces, and applications arise in 
adhesives, magnetic disks and membranes. They have 
therefore been the focus of many experimental and theo- 
retical studies @- When forced to cover a non-wettable 
substrate, a thin liquid film is unstable and will dewet 
this substrate. Four years ago, Reiter studied the dewet- 
ting of ultrathin (i.e. thinner than the coil size), al- 
most glassy polystyrene (PS) films deposited onto sili- 
con wafers coated with a polydimethylsiloxane (PDMS) 
monolayer Q. He found that a highly asymmetric rim, 
with an extremely steep side towards the interior of the 
hole and a much slower decay on the rear side, builds 
up progressively. In order to explain these deviations 
from the behavior of simple Newtonian liquids, theoret- 
ical models based on the shear-thinning properties of 
polymer films |5j has been proposed by Saulnier et al. 

and Shenoy et al. ML Very recently, Damman, Re- 
iter and collaborators [a compared the opening of 
cylindrical holes with the retraction of a strait contact 
line. These experiments revealed the existence of a highly 
asymmetric rim even in the case of the edge geometry, 
and a strong decrease of the dewetting velocity as the rim 
builds up (for both the radial and the edge geometries) 
[To|. One can show that none of these features can be 
explained by the existing theoretical models. In this let- 
ter, we therefore present a new theoretical approach that 
shows how the friction (due to slippage) of the liquid 
film onto the substrate, combined with the viscoelastic 
properties of the polymer film, can simply explain the 
above-mentioned experimental results of Damman and 
Reiter [10. 

Let us consider the dewetting of a straight contact line 
and assume that the thickness of the liquid film is smaller 
than the hydrodynamic extrapolation length b = rj/Q 
(where rj is the viscosity of the liquid, and £ the fric- 
tion coefficient of the liquid on the substrate) 0. We 



FIG. 1: Film geometry : h(r,t) is the profile of the film, ho 
is the initial height of the film, H(t) is the height of the front, 
L(t) is the dewetted distance, W(t) is the width of the rim, 
and v(x,t) is the velocity of the film. 



can thus use a simple plug-flow description and charac- 
terize the velocity field in the film, v(x,t), and the film 
profile, h(x,t), by two functions independent of the z- 
coordinate (see Fig. [IJ. The horizontal stress, a(x,t), is 
related to the strain rate, 7 = dv/dx, by a constitutive 
equation (the form of which depends on the type of fluid 
under consideration; see below). Neglecting inertia, local 
mechanical equilibrium between the friction forces (per 
unit surface) onto the substrate, ^v(x,t), and the bulk 
viscous forces gives: 



d(hcr) 

dx 



(1) 



Assuming the fluid to be incompressible, volume conser- 
vation leads to: 



dh 



dh 

dx 



dv 

dx 



(2) 



The last relation needed in order to solve the problem 
is the boundary condition at the edge of the film. The 
applied force (per unit of length) on the rim , \S\ (where 
S is the spreading parameter |lj assumed to be nega- 
tive), pushing the film away from the dry area, must be 
balanced by the the viscous force: 



\S\ 



-Ha(x = L) 



(3) 
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FIG. 2: Numerical calculation of the shape of a Newtonian 
film dewetting on a slippery substrate at different times. The 
different curves are separated by the same time interval. 

where H = H(t) is the front height, and L = L(t) is the 
dewetted distance (see Fig. 

For a Newtonian liquid, the above equations can eas- 
ily be solved at short times. Indeed, as long as h(x,t) 
remains of the same order as the initial thickness of the 
film, ho, equation Q - combined with the fact that for a 
Newtonian fluid a = 777 - leads to: 

, d 2 v ... 
C"** h o Wi (4) 

The velocity field is then given by v(x) = 
Vo exp (— ^^ ± ), where the distance A is given by 
\Jh a rj/Q = ^hob, and the velocity Vq by \S\/y/£rjh~o~. At 
short times, a Newtonian liquid deposited on a slippery 
substrate thus dewetts with a constant velocity Vq. 
This result was already obtained by Brochard-Wyart 
et al. 11] using energetic arguments, but the present 
mechanical point of view gives us additional informations 
about the film morphology. Indeed, Eqs. (J2J and @ 
give for the front height : H = ho + (\S\/rj) t. Since the 
velocity field decreases exponentially as one moves away 
from the front, the film profile exhibits at short times 
an asymmetric rim, with an exponential decrease of the 
thickness over the characteristic length A: 

h{x,t) = h + y* ex p( — 2Ty (5) 

This behavior is indeed the observed by Reiter on AFM 
images |4|. We have also completed our analysis by nu- 
merically solving Eqs. QJ, J2J and J3J; as shown on Fig. [21 
these numerical solutions confirm well our analytical pre- 
dictions. Our analysis also allows us to correct an as- 
sumption made by Brochard-Wyart et al. stating that 
the viscous dissipation should be negligible compared 
with dissipation due to friction [TH ]. Indeed, a simple 
calculus based on the above results shows that the two 
dissipations are approximately equal. Note that due to 
surface tension, the rim is in fact rounded over a distance 



S ~ H/6q (where 9q is the equilibrium contact angle). 
But the rim remains highly asymmetric if S <C A, that is 
to say as long as if A 7 or, equivalently, as long as 
L <C 9ob. For L > Oob, the friction of this cylindrical sec- 
tion on the substrate begins to be more important than 
the friction of the rest of the film, and simultaneously the 
Laplace pressure due to the curvature of the surface of the 
rim becomes stronger than the capillary pressure \S\/ H . 
Thus, once L > 80b, the "mature rim" regime described 
by Brochard-Wyart et al. ^HG2 (see also Damman et 
al. begins, and the rim becomes round and symmet- 
ric (with a width W ~ \fh~oL simply given by volume 
conservation). In this regime, the viscous dissipation is 
negligible compared with the dissipation due to friction, 
and, consequently, the dewetting velocity is proportional 
to 11]. Two important results arise from our anal- 
ysis of the dewetting of a Newtonian fluid. Firstly, the 
friction of the film onto the substrate gives rise to an 
asymmetric rim, since it dumps the velocity field in the 
film over a length A (which depends on the liquid viscos- 
ity). Secondly, the viscous dissipation is approximately 
equal to interfacial dissipation due to friction during the 
formation of the rim, while it is negligible in the "mature 
rim" regime. We can therefore anticipate that for a vis- 
coelastic fluid, the rheologic properties of the fluid will 
have no significant consequences on the dewetting veloc- 
ity in the "mature rim" regime, but will play a major 
role during the formation of the rim. 

Let us now consider in some details the dewetting of a 
viscoelastic film, assuming the following simplified con- 
stitutive equation |l3j: 

Go + (rjo + 771)0- = Grj! 7 + 770771 7 (6) 

where G is an elastic modulus (due to entanglements), 
?7o is a short time viscosity and 771 is the usual melt 
viscosity (771 3> 770). The time response of such a liq- 
uid can be divided into three regimes: (1) At short 
times, t < To = 770/G, the liquid behaves like a sim- 
ple Newtonian liquid with weak viscosity 770; (2) For 
To < t < n = rji/G, where t\ is the relaxation time of the 
liquid {i.e. the reptation time of the polymer chains) the 
liquid behaves like an elastic solid of elastic modulus G. 
(3) At long times (t > Ti), the liquid behaves like a very 
viscous Newtonian liquid of viscosity 771 . The above men- 
tioned time response of the liquid has direct consequences 
on the dewetting process. For times shorter than to, the 
viscoelastic liquid dewets like a simple liquid, with a con- 
stant velocity Vo — \S\/ y/^rjoho, and with the formation 
of an asymmetric rim of width Ao = yhorjo/C- At long 
times (t > n), the viscoelastic liquid also dewets like a 
simple liquid, with a constant velocity V\ = \S\/ y/QryJio, 
and with the formation of an asymmetric rim of width 
Ai = yho^i/C, ^ Ao. In between these two regimes, 
the viscoelastic behavior of the fluid will thus lead to 
a significant drop of the dewetting velocity (from Vo to 
Vi). More precisely, in this intermediate time regime 
(tq < t < Ti), the liquid behaves like an elastic solid and 
the height of the front increases very slowly with time: 
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FIG. 3: Numerical calculation of the reduced dewetting ve- 
locity V/Vo versus the reduced time t/ro for a viscoelastic film 
with ti = 1000 tq. The straight line represents (t/ro)^^ . 



H ~ h Q + \S\ (1 + t/n) /G ps h + \S\/G [lj. Volume 
conservation then imposes that the width of the rim, W, 
increases proportionally to the dewetted distance, L: 



W 



G_ 

W\ 



h L 



(7) 



Let us now assume that the bulk viscous dissipation is at 
most of the order of the dissipation due to friction. This 
assumption then allows us to determine the dynamic of 
the dewetting process from a simple energy balance (per 
unit of length) between the work done by the capillary 
force per unit of time and the dissipation due to friction: 



\S\V ~ (WV 2 



(8) 



(where, as anywhere else in this letter, numerical pref- 
actors of order unity are neglected). The above equa- 
tion, combined with Eq. Q), gives V(t) ~ Vq \J t / t for 
times t shorter than t\. Note that this t~i behav- 
ior might h ave be en guessed directly from the fact that 
Vi/Vq = \/To/ri. We have confirmed these results by 
numerically solving the equations of motion (see Fig- - 
The good agreement between our analytical and numer- 
ical results indicates that our assumption that the bulk 
viscous dissipation is smaller (or equal) to interfacial dis- 
sipation does hold for a viscoelastic liquid described by 
Eq. © • For t > t\, our model predicts a constant dewet- 
ting velocity V%. This prediction holds as long as the 
height of the front H is small compared to the width 
8qAi. Thereafter, as in the case of a Newtonian liquid, 
the " mature rim" regime is reached and the dewetting ve- 
locity then decreases like t~ 3 . We have thus shown that 
for a viscoelastic fluid in the time interval To < t < t\ , the 
width of the rim W increases proportionally to the dewet- 
ted distance L, and the dewetting velocity decreases like 
t~ 2 . The former prediction is in good agreement with the 
experimental observations of Damman et al. @ . While 
aT5 decrease of the dewetting velocity has indeed been 



occasionally observed by Damman [l5j , in most cases the 
measured dewetting velocity decreases like i _1 0, much 
faster than predicted. In order to explain this striking 
result, let us assume that at the beginning of the dewet- 
ting process the polymer films are not at equilibrium, 
and display residual stresses due to the spin-coating fab- 
rication process and fast evaporation of the solvent, as 
recently emphasized by Reiter and De Gennes [lfj. We 
shall now show how these residual stresses, assumed to be 
essentially horizontal and of initial amplitude <t , cause a 
high initial dewetting velocity, followed by a strong slow 
down. 

The various time regimes of the dewetting process, 
when partially driven by residual stresses, are similar 
to the ones already described when the process is only 
driven by capillary forces. At times shorter than To, 
the dewetting velocity is equal to Vq + <7oAo/?7o, where 
the second term denotes the contribution of the resid- 
ual stresses. In this short times regime, the residual 
stresses have no direct effects on the shape of the rim 
which thus keeps an exponential shape of characteristic 
width Ao. For t > tq, the height of the front is given 
by H = h + h cT /G + \S\(1 +t/r 1 )/G [HI, and the 
width of the rim, W, is simply given by volume conser- 
vation: W(H — ho) = hoL. Thus, as long as t <C ti, 
H is approximately constant and W increases propor- 
tionally to the dewetted distance. In order to obtain the 
dynamics of the dewetting process, the power (per unit 
of length) ho<jQ exp (— ~) V 01 delivered by the resid- 
ual stresses should be added to the l.h.s. of the energy 
balance Eq. © . The dewetting velocity is then given by 
(t>T ): 



V 



(1 



-)(1 



V2 



(9) 



where e = ho<Jo/\S\. Around t = To, the velocity de- 
creases like t^^, and thereafter decreases more sharply 
as the residual stresses relax in the film. For large enough 
residual stresses (e ^ 4), the dewetting velocity behaves 
like £ -1 around t = 2ti/3. Note that when the capil- 
lary forces are negligible (i.e. when £ > 1), the resid- 
ual stresses alone are able to induce the dewetting pro- 
cess and lead to a decrease of the dewetting velocity like 
exp (—t/n) (in the range n < t < t% In (e(l + e)). 

The above analytical results are in good agreement 
with numerical solutions (see Fig. 0}. Again, the simpli- 
fied energy balance resulting from our assumption that 
the bulk dissipation is smaller (or equal) to the interfa- 
cial dissipations gives very satisfying results. Residual 
stresses are thus a very good candidate to explain the 
experimental observations of Reiter and Damman since 
both morphological observations (rim width W propor- 
tional to the dewetted distance L), and dynamic mea- 
surements (variations with time of L and V) are in very 
good agreement with the theoretical predictions. Addi- 
tionally, it has been observed that the dewetting velocity 
of a circular hole, while much lower than the dewetting 
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FIG. 4: Numerical calculation of the reduced dewetting ve- 
locity, V/Vi, versus reduced time, t/ro, for a viscoelastic film 
(t~i = 100 to) with residual stresses (ao = G = 4\S\/ho). 
The straight lines represents (t/ro) _1//2 and (t/ro) -1 respec- 
tively. 

velocity of a strait contact line at the beginning of the 
dewetting process (due to radial deformations), system- 
atically joins it after some time 0, |tj . The presence of 
residual stresses can simply explain this experimental ob- 
servation. Indeed, as described above, the dewetting ve- 
locity is mainly controlled by residual stresses (for e > 1) 
which are evenly distributed throughout the film. Thus, 
when the size of a hole is large enough for the viscous 
dissipations due to radial deformations to be negligible 
compared with the friction onto the substrate, the hole 
becomes equivalent to a strait line, and both velocities 
become of the same order, even-though the dewetted dis- 
tances and the rim sizes arc different in both cases. 



In conclusion, we have shown that the friction (due 
to slippage) of the liquid film onto the substrate can ex- 
plain the building up of the asymmetric rim observed 
by Reiter || during the dewetting of thin PS films on 
a PDMS monolayer. We have also shown that the vis- 
coelastic properties of the PS are of great importance 
as they lead to a decrease of the dewetting velocity with 
time proportional to t~ 2 for times shorter than the repta- 
tion time of the polymer chains in the film. This decrease 
is made sharper by the presence of residual stresses. A 
sharp decrease of the dewetting velocity (V ~ t^ 1 ) as 
observed by Damman et al. 8] could thus be seen as an 
evidence of the presence of residual stresses in such vis- 
coelastic films. Note that these residual stresses should 
also influence the initial stage of the opening of cylindri- 
cal holes (a situation where the dissipation due to radial 
deformations dominates over the friction). The residual 
stresses might also play a role in the surface instabilities 
of the film and in the rate of hole formation ^(| . In this 
letter we did not talk about the shear thinning properties 
of PS films [HEl, but one can show - using an analysis 
similar to the one used in this letter for viscoelastic film 
- that a shear-thinning behavior leads to a decrease of 
the dewetting velocity weaker than t~ 2. Hence, in the 
absence of residual stresses, shear thinning alone cannot 
explain the obeservations of Reiter, Damman and colab- 
orators, even combined with viscoelastic properties. 



Acknowledgments 

We wish to acknowledge very interesting and fruitful 
discussions with Giinter Reiter and Pascal Damman. 



[1] P.-G. de Gennes, F. Brochard-Wyart and D. Quere, Cap- 
illarity and Wetting Phenomena: Drops, Bubbles, Pearls, 
Waves (Springer, 2003). 

[2] A. Oron, S. H. Davis, and S. G. Bankoff, Rev. Mod. Phys. 
69, 931 (1997) 

[3] G. Reiter, Science 282, 888 (1998) 

[4] G. Reiter, Phys. Rev. Lett. 87, 186101 (2001) 

[5] K. Dalnoki-Veress, B.G. Nickel, C. Roth, and JR. 
Dutcher, Phys. Rev. E 59, 2153 (1999) 

[6] F. Saulnier, E. Raphael, and P.-G. de Gennes, Phys. Rev. 
Lett. 88, 196101 (2002); F. Saulnier, E. Raphael, and P.- 
G. de Gennes, Phys. Rev. E 66, 061607-(1-12) (2002) 

[7] V. Shenoy and A. Sharma, Phys. Rev. Lett. 88, 236101 
(2002). For another theoretical approach, see S. Herming- 
haus, R. Seemann, and K. Jacobs, Euro. Phys. J. E 87, 
196101 (2001) 

[8] P. Damman, N. Baudelet, and G. Reiter, Phys. Rev. Lett. 

91, 216101 (2003) 
[9] G. Reiter, M. Sferrazza, and P. Damman, Euro. Phys. J. 

E 12, 133 (2003) 
[10] Note, however, that the geometry has drastic influence 

on the dynamics at early stage; see || 



[11] F. Brochard-Wyart, G. Debregeas, R. Fondecave, and P. 
Martin, Macromolecules 30, 1211 (1997) 

[12] C. Redon, J. B. Brzoska, and F. Brochard-Wyart, Macro- 
molecules 27, 468 (1994) 

[13] R. B. Bird, R. C. Armstrong, and O. Hassager, Dynamics 
of polymeric liquids (John Wiley & Sons, volume 1, 1977) 

[14] This last formula is valid assuming \S\ < hoG. More gen- 

erally, on has H = Ho + with % In (f^) = 

[15] P. Damman, private communication. 
[16] G. Reiter, and P.-G. de Gennes, Eur. Phys. J. E. 6, 25 

(2001) 

[17] This last formula is valid assuming \S\ + hoao < hoG. 
A more general formula can be obtained along the lines 

of [3. 

[18] The power delivered by the residual stresses is formally 
given by J h(x,t)a(x,t) A /(x,t)dx. Using the fact that 
a(x,t) vanishes at the front (x — L), and reaches it 
unperturbed value ao exp ( — — ) far away from the rim 
(x 3> W), one can obtain the simplified expression men- 
tioned in the text. 



